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Abstract 

Substantial progress has recently been reported in the determination of the Hilbert-Schmidt (HS) 
separability probabilities for two-qubit and qubit-qutrit (real, complex and quaternionic) systems. 
An important theoretical concept employed has been that of a separability function. It appears 
that if one could analogously obtain separability functions parameterized by the eigenvalues of 
the density matrices in question-rather than the diagonal entries, as originally used-comparable 
progress could be achieved in obtaining separability probabilities based on the broad, interesting 
class of monotone metrics (the Bures, being its most prominent [minimal] member). Though large- 
scale numerical estimations of such eigenvalue-parameterized functions have been undertaken, it 
seems desirable also to study them in lower-dimensional specialized scenarios in which they can be 
exactly obtained. In this regard, we employ an Euler-angle parameterization of 50(4) derived by S. 
Cacciatori (reported in an Appendix)-in the manner of the 5'C/(4)-density matrix parameterization 
of Tilma, Byrd and Sudarshan. We are, thus, able to find simple exact separability (inverse-sine- 
like) functions for two real two-qubit (rebit) systems, both having three free eigenvalues and one 
free Euler angle. We also employ the important Verstraete-Audenaert-de Moor bound to obtain 
exact HS probabilities that a generic two-qubit state is absolutely separable (that is, can not 
be entangled by unitary transformations). In this regard, we make copious use of trigonometric 
identities involving the tetrahedral dihedral angle (p = cos^^ (^). 
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I. INTRODUCTION 

Zyczkowski, Horodecki, Sanpera and Lewenstein (ZHSL) P] were the first, it is clear, 
to pose the interesting question of determining the probability that a generic two-qubit or 
qubit-qutrit state is separable [1] . The present author has further pursued this issue p]4TT] . 
In particular, he has addressed it when the measure placed over the two-qubit or qubit- 
qutrit states is taken to be the volume element of certain metrics of interest that have been 
attached to such quantum systems (cf. [121 sec. 5]). (This parallels the use of the volume 
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element ["Jeffreys' prior] of the Fisher information metric in classical Bayesian analyses 
[131 [II].) The (non-monotone [12] ) Hilbert-Schmidt [TE] and (monotone [I3|) Bures [IH] can 
be considered as prototypical examples of such quantum metrics [19] chap. 14]. 

One useful (dimension-reducing) device that has recently been developed, in this regard, 
is the concept of a separability function. (The dimension-reduction stems from integrating 
over-that is, eliminating- certain [off-diagonal] parameters.) The integral of the product of 
the separability function and the corresponding jacobian function yields the desired separa- 
bility probability ^ eqs. (8), (9)]. In the initial such studies [HI [ID], the separability functions 
were taken to be functions of (ratios of products of) the diagonal entries of the 4 x 4 or 
6x6 density matrix (p). In the later paper [llj, sec. Ill], it was also emphasized that it 
would be desirable to find analogous separability functions parameterized, alternatively, in 
terms of the eigenvalues of p. The motivation for this is quite straightforward. "Diagonal- 
entry-parameterized separability functions" (DSFs) have proved quite useful in studying 
Hilbert-Schmidt separability probabilities. However, the formulas [IH] for the large, inter- 
esting class of monotone metrics (Bures /minimal monotone, Kubo-Mori, Wigner-Yanase,...) 
are expressed in terms of the eigenvalues, and not the diagonal entries of p. Until such 
"eigenvalue-parameterized separability functions" (ESFs) are obtained, it does not appear 
that as much progress can be achieved as has been reported for the Hilbert-Schmidt metric 
employing DSFs [9l-[TT]. (Let us take note here-although we are not aware of any imme- 
diate relevance for the problems at hand-of the Schur-Horn Theorem, which asserts that 
the increasingly-ordered vector of diagonal entries of an Hermitian matrix majorizes the 
increasingly-ordered vector of its eigenvalues [201 chap. 4] (cf. [21]).) In particular, we 
would hope to be able to test the conjecture [numerically suggested] that the two-qubit 
Bures separability probability is ^680(v^-i) ^ 0.0733389 [6]. 

In [H] we also undertook large-scale numerical (quasi-Monte Carlo) analyses in order to 
estimate the ESF for the 15-dimensional convex set of two-qubit states. We have continued 
this series of analyses (also now for the 9-dimensional real two-qubit states). However, at 
this stage we have not yet been able to discern the exact form such a (trivariate) function 
putatively takes. In light of such conceptual challenges, it appears that one possibly effective 
strategy might be to find exact formulas for ESFs in lower-dimensional contexts, where the 
needed computations can, in fact, be realized. (This type of "lower-dimensional" strategy 
proved to be of substantial suggestive, intuition-enhancing value in the analyses of DSFs. 
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One remarkable feature found was that the number of variables naively expected [that is, 
n — 1, when p is n x n] to be needed in DSFs can be reduced by substituting new variables 
that are ratios of products of diagonal entries-so the individual independent diagonal entries 
need not be further utilized. One of our principal goals is to determine whether or not 
similar reductive structures are available in terms of ESFs.) 

Tilma, Byrd and Sudarshan have devised an {SU{A)) Euler-angle parameterization of 
the (complex) two-qubit states [22]. To simplify our initial (lower-dimensional) analyses, 
we have corresponded with Sergio Cacciatori, who kindly developed a comparable (5*0(4)) 
parameterization for the 9-dimensional convex set of two-qubit real density matrices. (This 
derivation is presented in Appendix I.) In this 5'0(4)-density-matrix-parameterization, there 
are six Euler angles and three (independent) eigenvalues, Ai, A2, A3, with A4 = 1 — Ai — A2 — A3 
(cf. m eq. (9)]). 



II. DETERMINATION OF ESFS 



A. Example 1 

To this point in time, we have only been able to obtain exact ESFs when no fewer than 
five of the six Euler angles are held fixed. We now present our first such example, allowing 
the Euler angle xi to be the free one, and setting (using, for initial simplicity, the midpoints 



P 



of the indicated variable ranges (49)) X4 = Xq = ti and X2 = x^ = x^ = |. We will, thus, to 
begin with, be studying density matrices of the form, 

/ Ai 

cos^ (xi) A2 - sin^ (xi) (Ai + A2 + A3 - 1) sin (2xi) (Ai + 2A2 + A3 - 1) 
-|sin(2xi)(Ai + 2A2 + A3- 1) sin^ (xi) A2 - cos^ (xi) (Ai + A2 + A3 - 1) 

Vo 

(1) 

By the Peres-Horodecki positive-partial-transpose (PPT) condition [211125], P is separable 
if and only if (an irrelevant nonnegative factor being omitted from the PPT determinant) 

A1A3 - ^ sin2(2xi) (Ai + 2A2 + A3 - 1) ' > 0. (2) 



Since we have set X2 = X3 = X5 = ^, the Haar measure (48) simply reduces to unity. 



Integrating this measure over the interval Xi G [0,27r], while enforcing the separability 



condition (|2])-as well as requiring the facilitating eigenvalue-ordering Ai > A2 > A3 > A4-we 
obtain the desired ESF{Xi, A2, A3). 

This ESF is equal to unity under the constraints (recall that we set A4 = 1 — Ai — A2 — A3) 

Ai + 2A2 > 1 A 2A2 + A3 + 2v/Ai - 2A1A2 > Ai + 1 (3) 

and 

Ai + 2A2 = 1 V 2A2 + A3 + 2^\i - 2A1A2 > Ai + 1. (4) 

These constraints, thus, define the domain of eigenvalues for which all the possible density 
matrices of the form ([T])-independently of the particular value of xi E [0, 7r]-are separable. 
Now, nontrivially, outside this domain of total separability ([s]), (|4]), we have 

E5F|,,}(Ai,A2,A3) = ^^^^, (5) 

where 

^-T^ (6) 

A2 — A4 

and 

A1 + 2A2 > 1 A2A2 + A3 + 2v/a7^2A^< Ai + 1. (7) 

(The set of constraints ([T]), together with the imposed nonascending order of the eigenvalues, 
ensure that the argument of the inverse sine function, U, is confined-sensibly-to the interval 
[0,1].) The argument U can be seen to be intriguingly analgous to the important (abso- 
lute) separability bound-originally suggested by Ishizaka and Hiroshima [26]-of Verstraete, 
Audenaert and De Moor [27] [2H1 eq. (3)], 

VADiX,, A2, A3) = Ai - A3 - < 0, (Ai > A2 > A3 > A4), (8) 

or, equivalently, 

C7.^>1. (9) 

Ai — A3 

(The associated greatest possible value of Ai is, then, ^; of A2, ^-'g^; of A3, |; and of A4, |.) 
The inequality ^ improves upon the ZHSL purity (inverse participation ratio) bound [T] 
(Figs, g and § (cf. ([MEI])), 

Tr(p2) = \i + \l + Xl + \l<l. (10) 

It appeared to us that the ESF for the full 9-dimensional real and/or 15-dimensional 
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FIG. 1: Three-dimensional contour plot of F^D(Ai, A2, A3), given by eq. ([s]). Axes are of length 
1. 

complex convex sets of two-qubit states might be of the form, 

TT 

for VAD{Xi, A2, A3) > 0, and simply unity for VAD{Xi, A2, A3) < 0. 
B. Example 2 

For our second example with one free Euler angle and three free eigenvalues, we set 
Xi = X4 = xq = n and X2 = X3 = ^, allowing X5 to be the free Euler angle now. (Scenarios 
with only one of 0:2, 0:3, 0:4 or xq allowed to be free are trivially completely separable, and 
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FIG. 2: The absolutely separable boundary surface VAD{Xi, A2, A3) = 0. Axes are of length ^. 

thus do not merit attention.) This results in density matrices of the form (again, taking 
A4 = 1 — Ai — A2 — A3), 

^ A3 cos^(x5) + sin^(a;5)Ai cos(x5) sin(x5)A3 — cos(x5) sin(a;5)Ai 

A2 

-Ai - A2 - A3 + 1 

\ cos(x5) sin(x5)A3 — cos(x5) sin(x5)Ai Ai cos^(x5) + sin^(x5)A3 

(12) 

The Peres-Horodecki separability condition can now be expressed as 

- ^ sin^ (2x5) (Ai - A3) ' - A2 (Ai + A2 + A3 - 1) > 0. (13) 



In contrast to our first example, however, the (reduced) Haar measure (48) is now non 



uniform, being sin 0:5. We, thus, sought to integrate this measure over the range G [0, tt], 
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FIG. 3: The surface corresponding to non- absolutely separable states with VAD{Xi, X2, \3) = \. 



Axes are of length ^ 



3- 

3- 



while enforcing the condition (13), to obtain the corresponding ESF. However, Mathematica 
did not yield a solution after what we judged to be a reasonable amount of computer time 
to expend. 

Therefore, we made the transformation X5 = cos~^ {—y^)-, which leads to a uniform 
measure for G [—1,1]. (There were some initial concerns expressed, in regard to this 
transformation, by S. Cacciatori. He had written: "The transformation you consider is not 
bijective on the whole range of the parameters. Indeed, the uniformization can be done 
only locally. This is because the measure on a compact manifold cannot be an exact form 
(otherwise, using Stokes Theorem, one finds zero volume for the manifold)" (cf. j32l p. 
4394]). However, subsequent analysis revealed that nothing fallacious arises in this manner, 
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FIG. 4: The blue points correspond to values of the (unsorted) eigenvalues for which the VAD 
bound on absolute separability ([s]) is equivalent to the ZHSL purity bound, + + A| + A| < |. 
The red points are those which satisfy the VAD bound, and for which the purity is, then, the 
greatest possible, that is, | > |. 

since the results are equivalent to twice those obtained by integrating over (the half-range) 
1/5 G [0, 1].) Now, the region of total separability was of the form, 

Ai < ^ A Ai > A ((2Ai + As > 1 A (M2 V M3)) V (M4 A M5)) , (14) 

where 

M2 = Ai + 2^/X2 - 2A1A2 = 2A2 + A3, M4 = 2A2 + A3 + 2^/X2-2X,X2 > Ai, (15) 
M3 = 2A2 + A3 + 2v/A2 -2A1A2 > Ai A Ai + 2^/X^^2XJ^2 > 2X2 + A3 (16) 

and 

M5 = (2Ai + A2 = 1 A Ai + A2 + A3 < 1) V (|A2 > A Ai + 2v^A2 - 2A1A2 > 2A2 + A3 A 2Ai + A2 < 1 

(17) 
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FIG. 5: Fig. [4] embedded in the simplex of possible values of three (unordered) eigenvalues 
Also, non-trivially, the ESF took the (root) form 

ESF{^^y{U) = Root [4#1^ - 4#l2 + U^L,3] - Root [4#1^ - 4#1^ + U^k,A] + 1 (18) 
(Root[/, k] represents the k-th root of the polynomial equation f[x] = 0) for 
2Ai + A2 > 1 A < Ai < ^ A Ai - 2A2 - A3 + 2^X2 - 2A1A2 < A Ai + A2 + A3 < 1 (19) 



or 



< Ai < 1 A Ai + A2 + A3 < 1 A A2 > 0. 



(20) 



In fact, (18) has an equivalent radical form. 



V2 



V2 



(21) 



The equivalence can be seen from a joint plot, in which ( 18 ) amd (21 ) fully coincide. However, 



the Mathematica command "ToRadicals" did not produce (21), but gave the sign of the 



second of the three addends as plus rather than minus. ("If Root objects in expr contain 
parameters, ToRadicals [expr] may yield a result that is not equal to expr for all values of the 
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parameters".) This was clearly an erroneous ("hyperseparable" -greater than 1) result. We 
initially had thought it might be due to an inappropriate uniformization, = cos"^ ("Z/s)- 



(In Appendix II, we also give a derivation of S. Cacciatori of (21) "by hand".) Therefore, 



we had also investigated alternative approaches to obtaining the required ESF. 

A1-A3 



We had reasoned that since the transformed VAD variable U = V^^^'* seemed to play a 



vital role, we performed the change-of-variables. 



|a2 ^ ^ (-Ai - A3 + v/(Ai + A3-l)2-F(Ai-A3)2 + 1) I , (22) 

in our constrained integration analyses, where V = U"^. Now, we again imposed the 
nonascending-ordering requirement on the eigenvalues (and their transformed equivalents). 



and the Peres-Horodecki condition (13) (again discarding irrelevant nonnegative factors) 
became 

2V + cos (4x5) - 1 > (23) 
(conveniently being free of the individual A's). Integrating the reduced Haar measure sinxs 



over X5 G [0,7r], while enforcing (23), we obtained the result (Fig. |6]), 

1 V>1 



1 - ^ cos {I sm-\l - 2V)) - ^1 + ^ sin (1 sm-\l - 2V)) +1 < V < 1. 

(24) 

We can represent this ESF more succinctly still, by using the variable W = 1 ~ 2V. Then, 
we have 

|l W<-1 
ESF{,,y{W) = { . (25) 

I 2sin2 (|cos-i(H^)) +sin {lcos-\W)) -1 < < 1 
III. ABSOLUTE SEPARABILITY ANALYSES 



A. ZHSL uniform simplex measure 

Verstraete, Audenaert and De Moor remarked that the use of their result ([s]) yields a 
better lower bound (assuming a uniform distribution on the simplex of eigenvalues) for the 
volume of separable states relative to the volume of all states: 0.3270 (as opposed to the 
ZHSL bound of 0.3023 based on the purity [H eq. (35)]) |27l p. 6]. We are now able to 
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FIG. 6: The eigenvalue-parameterized separability function given by eq. (24) — along with the 

2sm-l(v^) 



slightly subordinate function, (11) 



-, y G [0, 1] — for the real two-qubit scenario with the 



Euler angles xi = 7r,X2 = f,iC3 = |, X4 = vr, xe = tt and X5 free. Here, V = U'^ = p^^r^p") with 
Ai > A2 > A3 > A4 



give exact formulas for these two bounds (certainly the second being new, and apparently 
somewhat challenging to derive), namely 



ZHSL 



IT 



lower— bound 



0.30229989 



(26) 



and (Figs. Ill |2] and [2 



6v/2 + 3v/2tan- ( '-^^ 



VADi,^,r-bound = ^ I 8 - 6V2 + 3V2tan-^ I ^^^^ ) ) ~ 0.32723006. (27) 



The Verstraete-Audenaert-de Moor bound (|8j) is based on the entanglement of formation. 
They also present another bound [271 P- 2] 

VAD{Xi, A2, A3) = V(Ai - A3)2 + (A2 - A4)2 - A2 - A4 < 0, (Ai > A2 > A3 > A4), (28) 

based on the negativity. (There is also a more complicated bound based on the relative 
entropy of entanglement. Since it involves logarithms, it is more difficult for us to analyze.) 
From this we obtained. 



VAD, 



owei — bound 



12 - 8v^ + 37r + 12 cot-i ( ~'''™^ ) 
4^2 



0.32723006, 



(29) 



fully equivalent, it would appear, to (27). (Perhaps it is clear that these two bounds (pi) and 



(28) must be equivalent, though no explicit mention is made of this, it seems in |27j nor the 
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more recent [28]. Mathematica readily confirmed for us tliat tliere is no nontrivial domain 
in wfiicli one of tfie bounds is greater tfian and the otlier one, less than 0.) Let us observe 



further, in regard to (29) that 6561 = 3 . Apparently, there is a right triangle with opposite 
legs of lengths 5983 and 1904a/2 and hypothenuse 6561 = 3®, that is of relevance in some way 
to the issue of the absolute separability of two-qubit states. We have consulted the Integer 



Superseeker website of N. J. A. Sloane (http://www.research.att.com/~njas/sequences) in 
regard to this issue, and found sequence A025172, entitled " Let = arccos (|), the dihedral 
angle of the regular tetrahedron. Then cos (n0) = In our application, n = 8 and 

a{n) = —5983. The accompanying comment is that the sequence is "Used when showing 
that the regular simplex is not 'scisssors-dissectable' to a cube, thus answering Hilbert's 
third problem." Applying these facts to the problem at hand, we are able to obtain the 
simplified results 

VADion^er-bound = V ADio^er-bound = ^ (8 - 6^2 + 3^2 (tt - 8 csc-i(3))) ^ 0.32723006. 

(30) 

We have obtained an area-to-volume ratio of (cf. Fig. |2]) 

12(-l + 3y2cot-^ (2v/2)) 12(-l + 3x/2cot-^(2y2)) 

KArea/Vol = ^ ^ = ^ ^ ^ ~ 4.U0U4i0. 

Area/Voi g _ + (vT - 8 CSC"! (3)) 8 - 6v^ + 3 V^VT - 24v^ cot'^ (2 v^) 

(31) 



B. HS measure on two-qubit real density matrices 

Further, if we employ the Hilbert-Schmidt measure for two-rebit density matrices [T6l 
eq. (7.5)] on the simplex of eigenvalues, we can obtain (V. Jovovic assisted with several 
trigonometric simplifications involving the dihedral angle, cos^^(|) an exact lower bound 
(much weaker than the conjectured actual value of ~ 0.470588 (THl sec. 9.1]) on the 
Hilbert-Schmidt two-qubit separability probability. This is 

^^^Ht-w.-.oW = - 0.0348338. (32) 



Also (cf. 



where 



RHS-real ^ 1 f ^ N _ 12.489976122, (33) 

Area/voi i287(-6928 + 22057r) ^ 1^2; , k j 



34087768^2 - 2478673447r - 1292769261 cos"^ Q 
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and 

02 = 3226925088 cot^^ (^2^ + 3760128 tan^^ (272) + ^50082810 tan^^ (^^2 

The corresponding ratio based on the Bures metric is modestly larger (as will also be the 
case in the corresponding complex and quaternionic scenarios), that is, 14.582. Multiplying 



(33) by the radius-:j^-of the maximal ball inscribed inside the simplex of eigenvalues, we 



obtain a dimensionless ratio, 7 ~ 3.60555 [33i eq. (1)] (cf. [51]). 
C. HS measure on two-qubit complex density matrices 

We have also obtained an exact expression for the HS absolute separability probability 
of generic (complex) two-qubit states (using the indicated measure [Ml eq. (3.11)]), 



HS— lower— bound gr 

where 



^^^H5-tL.-.oW = 32614907904 ^ 0.0036582630543035 (34) 



V^i = 1959684729929728 - 1601255307608064^2 - 1529087492782080V27r 

and 

= 45247615492565918250^2 cot"^ (V2) - 22619730179635540245^2 sec^^ (3). 



(The conjectured HS [absolute and non-absolute] separability probability of generic two- 
qubit complex states is ~ 0.242424 ^TO], sec. 9.2].) Further, we have 

RTm''^ = -3840^^^ ^ 20.9648519, (35) 

■73 + T4 

with 

n = -5358569267936 + 33756573946095^271 - 270052591568760^2 cot"^ (Vt) , 
T2 = 11149704525960^2 cot^i (2V2) + 270052591568760^2 cot"^ (s + V2^ 



Tg = -1959684729929728 + 1601255307608064v^ + 1529087492782080^271, 

and 

n = -45247615492565918250^2 cot-i (^2^ + 22619730179635540245^2 sec"^ (3). 
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The corresponding dimensionless ratio is 7 ~ 6.05203, while the Bures analogue of (35) is 
23.7826. 

Employing the Bures measure [IHl eq. (3.18)] for generic two-qubit states on the simplex 
of eigenvalues, we have, numerically, a lower bound-based on the absolutely separable states- 
on the two-qubit Bures separability probability of 0.000161792, considerably smaller than 
the conjectured value of ^^soy-i) _ 0.0733389 [6J. 

D. HS measure on two-qubit quaternionic density matrices 

The HS absolute separability probability of the quaternionic two-qubit states is 

VAD%t-,^„_,„^„. = - 816946343106356485029888 ^"-^- " •'■•'•'•'039870347068. (36) 

(the conjectured absolute and non-absolute separability probability being g^^^QY25 ~ 
0.0773765 dH eq. (15)]), where 

Ci = -216449750678398795533760757497856 + 176860737736399592490919645937664^2, 

C2 = 279292548969739228073088142369304501839785^271 
(3 = -558572941247617043110461841280869072896000^2 cot"^ (v^^ 

C4 = 23637916932187025487103667523337320^2 cot-^ (^2\/2^ 
Cs = -16178155879591789043088455851252390200^2 cot"^ (s + v^J 
= -558589165778586158484606527963549721006600v^tan-i fv^^ 
Additionally, 

13 

jlHS-quat ^ ^. . . _ 37.9283799507, (37) 

Area/voi 3606947894919 168 V"""* ^ ^ ' ^ ^ 



where ^Z"?""* is given by ( 36 ) and 

V-i = -18147776040854148031593056 - 4720063928074960763823525^271, 

and 

^2 = -37760511424599686110588200y2tan~^ (9 - 7^2^ . 
The corresponding dimensionless ratio is 7 ~ 10.948980, and the Bures area-to-volume ratio 



counterpart to (37) is-slightly higher as in the real and complex cases-42.115. 
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1. Two-qubit Clement- Rag gio spectral condition 



Clement and Raggio have given a simple (linear) spectral condtion that ensures separa- 
bility of two-qubit states [301 Thm. 1], 



3Ai + V2X2 + (3 - V2)\3 < 2, Ai > A2 > A3. (38) 

Based on this, we readily found-using the ZHSL uniform measure on the simplex of 
eigenvalues-that a lower bound on the separability probability is ~ 0.235702, and the 
area-to- volume ratio is 6. 

Now if we employed the HS measure for the real density matrices, along with the spectral 



condition (|38|), the lower bound is ^^^tjn/^ ~ 0.0120065 and the area-to- volume ratio, 18. 



The corresponding two figures in the HS complex measure case are 0.00060239769 and 30, 
and, in the quaternionic case, 1.502473896 ■ 10~^ and 54. 

Clement and Raggio also gave a spectral separability condition applicable to any finite- 



dimensional system, but "much weaker" than (38) in the specific two-qubit case |30l Thm. 
2]. In the two-qubit case, it takes the form, 

3A3 + 3A4 > 1, As > A4. (39) 

Based on this constraint, the lower bound on the separability probability using the ZHSL 
uniform measure is | ~ 0.11111, while the area-to-volume ratio is 6. The associated separa- 
bility probability based on the HS measure on the real density matrices is ~ 0.00106691 
and area-to-volume ratio, 18. The corresponding HS complex measure values are ^^^^^qq^ ~ 
9.96592 ■ 10"® and 30, while the two quaternionic values are 2541865828329 ^ 8.59605 ■ 10~^° 
and 54. 

We see that the four area-to-volume ratios given are identical using either of the two 



bounds (38) and (39), that is, 6, 18, 30 and 54.. 



E. Qubit- Qutrit analyses 

The counterpart of the VAD bound (|8| for qubit- qutrit states was obtained by Hildebrand 
[281 eq. (4)] (Fig. g 



^(Ai,A2,A3,A4,A5) = Ai - A5 -20^ < 0, (Ai > A2 > A3 > A4 > A5 > Ag). (40) 
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FIG. 7: Three-dimensional projection of the sets of six (unsorted) eigenvalues, \i,= 1, ... ,6, for 
which the Hildebrand bound ( [40| ) on absolute separability of qubit-qutrit systems is equivalent to 
the ZHSL purity bound SfA^ < i (cf. Figs, g and |5|). 

We have also sought to obtain exact lower bounds on separability probabilities here, us- 



ing (40). Discouragingly, however, it seemed, after some initial analyses that computer 
memory demands may be too great to make any significant analytical progress on the full 
five-dimensional problem. (Nor have we been able to determine-as we have in the qubit- 
qubit case (after (|9|)-the maximum values that the Aj's can attain for absolutely separable 
states.) However, we then sought to reduce/specialize the problem to more computationally 
manageable forms. 

Firstly, we found that if we set Ai = |, then the absolute separability probability em- 
ploying the uniform measure on the resulting four-dimensional simplex is 0.00976679. (The 
corresponding value for the full five-dimensional simplex was reported by ZHSL as 0.056 [H 
eq. (35)].) However, this probability jumps dramatically to 0.733736, if we reduce Ai from 
|to|. 
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1. Clement- Rag gio spectral conditions 



We did succeed, however, in this qubit-qutrit case in terms of the "simple spectral con- 
dition" of Clement and Raggio [501 Thm. 2], which in our case takes the form, 

SAe + 5A5 > 1, A5 > Ae. (41) 

The resultant bound on the ZHSL separability probability (based on the uniform measure 
over the simplex) is, then, quite elegantly, ^ = 2~^, and the associated area-to- volume ratio 
is simply 15. Additionally, we found that the lower bound for the separability probability 
based on the HS measure for real density matrices is 

18989 



214748364800000 
Again, we have a simple area-to-volume ratio, 60. 



8.842442 • 10"'\ (42) 



IV. DISCUSSION 

We have attempted, so far without success, to obtain exact ESFs with more than one 
free Euler angle. In particular, we have tried to "combine" the two one-free-Euler-angle 
scenarios analyzed above, by letting both xi and X5 be free. The associated Peres-Horodecki 
separability criterion (cf. (|2|, (13)) is then. 



(A3 cos^ (X5) + sin^ (xs) Ai) (Ai cos^ (xg) + sin^ (xg) A3) sin^ (2xi) (Ai + 2A2 + A3 - 1) ^ > 0. 

(43) 

The best we were able to achieve, having first integrated over X5 G [0, tt], was a finding that 
the ESF for this two-Euler-angle scenario (again assuming Ai > A2 > A3 > A4) is expressible 
as 

/•2^ I 1 K>1 
E5F|,„,,}(Ai,A2,A3)= / { dx, (44) 

cos (I cos^^{k)^ — sin (| cos^^{k)) —1<k<1 







where 



cos (4xi) (Ai + 2A2 + A3 - 1) ^ + 4A2 + 2A3 - 4A2 (A2 + A3) + Ai (-4A2 + 4A3 + 2) - 1 

(Ai-A3)2 

(45) 
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Alternatively, we were able to reduce this two-free-Euler-angle problem to 

f-^ j 1 rj > 1 

ESF{^^^^^y{Xi,X2,X3) = / sin 0:5 < _ dx^ (46) 

Jo ['-^ + 1 -1<V<1, 

where 

^ - cos (43:5) (Ai - A3) ' + 4A2 + 2A3 - 4A2 (A2 + A3) + Ai (-4A2 + 4A3 + 2) - 1 
^ (Ai + 2A2 + A3-1)2 • ^ ^ 

On the other hand, we have not yet found an exact ESF with even just one free Euler 
angle using the Tilma-Byrd-Sudarshan SU{A) parameterization ^22J of the 15-dimensional 
convex set of (in general, complex) 4x4 density matrices. Further, we have been able 
to convince ourselves-somewhat disappointingly-that the ESF for the real 4x4 density 



matrices can not simply be a (univariate) function of V (or f/)-such as (11). We were able 
to reach this conclusion by finding distinct sets of eigenvalues (Ai, A2, A3, A4), which yielded 
the same identical value of V (we used V = but were capable of giving opposite signs for 
the determinant of the partial transpose, when sets of six (randomly-chosen) Euler angles 
were held fixed. 

The analyses we have presented above have considerable similarities in purpose with a 
notable study of Batle, Casas and A. and A. R. Plastino, "On the entanglement properties 
of two-rebit systems" [23] (cf. [35] )• Some differences between our work and theirs are 
that we have employed the Verstraete-Audenaert-de Moor entanglement measure ^ and 
its transforms {U and V) rather than the purity (which gives a smaller domain of absolute 
separability), and we have striven to obtain exact analytical results, rather than numerical 
ones (cf. [m sec. 3]). (Let us note, in passing-although we have not adopted this viewpoint 
here- "that there is qualitative difference between the separability problems for real and 
complex matrices. In fact, the set of separable density matrices has the same dimension 
as the full set of density matrices in the complex case, but has lower dimension in the real 
case" [Ml p. 713].) 

We have been principally concerned here (sec. |ll]) with the generation (by integration 
of the Haar measure over the Euler angles) of metric-independent eigenvalue-parameterized 
separability functions (ESFs). Such ESFs could be applied, in conjunction with the formulas 
for arbitrary metrics (e. g. ^ eq. (4.1)], ^Si eq. (3.18)] [HI eqs. (14.34), (14.45)]) to yield 
separability probabilities specific to the metrics in question. 
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We have also reported (sec. Ill) advances in the study of absolute separability, imple- 
menting the Verstraete-Audenaert-de Moor bound ([s]) . A substantial analytical challenge in 
regard to this is the determination of the volumes and bounding areas (and their ratios) of 
the absolutely separable two-qubit states in terms of the Bures (minimal monotone) metric. 
(We have given nnmenca/ values of these ratios above.) We have derived several remarkable 
exact trigonometric formulas pertaining to absolute separability. (V. Jovovic has assisted 
us in this task, simplifying Mathematica output by using identities involving the tetrahe- 
dral dihedral angle = cos^^ (I)-) Nevertheless, we can certainly not ensure that further 
(conceivably substantial) simplifications exist. 
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V. APPENDIX I (S. CACCIATORI) 



The Lie algebra so(4) of 5'0(4) consists of the 4x4 antisymmetric real matrices. A basis 
is given by the following matrices: 



T2^ 



( ° ^ 





o\ 


-1 

















y 





0/ 


/ 


1 


o\ 











-1 








^ 







/o 





o\ 





1 





-1 








^0 





oy 


/ 





l\ 




















[-1 





oy 


/o 





o\ 








1 











[o -1 





oy 


/o 





o\ 

















1 


\0 - 


-1 


oj 



. The maximal proper subgroup of SO (A) is SO (3) and can be thought of as the isotropy 
subgroup of SO{4:). This means that if we consider 5'0(4) as the group of rotations of R^, 
then this group acts transitively on the unit sphere (translating the north pole everywhere 
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on the sphere), but any point is left fixed by the subgroup SO{3). Thus we will have 

SO{4)/SO{3) ~ S\ 



This is the main observation needed to determine the range of the parameters. 
Now let us look at the construction of the Euler parametrisation of the group. First, one 
searches for a maximal subgroup, in our case H — SO {3). (In general, this might be 
not unique: for example for Eq we have chosen F4 as a maximal subgroup, but another 
possible choice could have been C/(l) x Spin{10)). We can suppose that we know the 
Euler parametrisation for the subgroup. Otherwise, we can proceed inductively, choosing a 
subgroup for H and so on. 

In our case for H, we can choose the 5*0(3) subgroup generated by Tj, i = 1, 2, 3. The Euler 
parametrisation ol H — SO {3) is well known, 

h[x,y,z] = e^^^e^^^e^^^ . 

To determine the range of the parameters we first note that 5*0(3) is not simply connected, 
but its universal covering is SU{2) and 5*0(3) ~ SU{2)/{I,—I} (/ being the identity). 
Also, it is well known that SU (2) ~ 5*^ as a manifold. With this in mind, let us compute 
the invariant metric (and measure) on H. This can be done, noting that the Lie algebra is 
isomorphic to the tangent space at the identity e of the group. But our algebra is provided 
by a natural scalar product, the trace product 

( , ) : so(4) X so(4) — > R, 
{A,B) I — y -^Trace(AS) . 

It is symmetric and positive definite, satisfies 

and induces a Euclidean scalar product on Tgi?, with {Ti^Tj) — Sij. Thus, it is easy to 
compute the metric on H as the metric induced by the scalar product on TeH. Indeed, if 
h[x, y, z] is a generic point, then dh is cotangent at the point h and h~^dh is cotangent at e, 
so that the metric tensor can be defined as 

dsjj = {h-^dh,h-^dh) = - ^Treice{h-^dhh-^dh) . 

23 



Prom this expression, we immediately see that ds^ is both left and right invariant. A direct 
computation gives 

ds"^ — dx^ + dy^ + dz^ + 2 cos{y)dxdz. 

Now we need to compare this with the metric of a sphere S^. This can be done, describing 
as the unit sphere in C^, 

= {{z,w) e C'^\zz + ww = 1}. 

Then, the point p of can be parametrised as 

p= (cosAe*('^+'^),sinAe'(^-^)) 

covered by A G [0,7r/2], e [0,7r], ijj e [0,27r]. The Euclidean metric is 

ds^ = dzdz + dwdw = d\^ + + + 2 cos(2A)d0#- 

Comparing this with the invariant metric, we see that it is the metric of with y — 2A, 
X — 2(f) and z — 2ip (note that x and z can be symmetrically interchanged). Now if we 
impose on x,y,z the ranges for the sphere, that is 

xe[0,27r], ye[0,7r], ze[OA^], 

we will cover SO{3) exactly two times. This is because SO{3) ~ S^/±I as noted before. 
However, we can easily check where duplication takes place. Indeed, 

/l o\ 

™ cos 2; sin 2; 
—sin 2; cos 2; 
^0 1/ 

and takes all possible values twice when z varies in [0, Air]. The correct range is then 

xe[0,27r], Z/e[0,7r], ze[0,2TT]. 

Now we are ready to construct the group G = SO (A). Any element of G can be written in 
the form 

e^^'+y^'+'^'h[x^,X5,Xe]. 
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This is shown in an appendix of our paper on F4 [37]. Here z) varies in M^. (We are 
not concerned at this level with the possibility of covering the group many times. The only 
important point is that this parametrisation is surjective). Now, we note that H acts on T^, 
as the rotations of 5*0(3) act on the points of the unit sphere S"^. In particular, from 

[T2, T4] = -Te, [T2, n] = T4, [T3, T4] = 0, [T3, n] = T5, [T3, T5] = -Tg 

we find 

e^'^^T^e-'"'^^ = C0SXT4 - sinxTe , 

^yTs ^xT2rp^^-xT2 ^-yTs _ (.Qg 3, _ sin x COS j/Tg — sinxsin^/Ts. 
Because any vector in can be written as 

X = r(cosx, — sin X sin?/, — sinx, cosy), 

and because 
we can write 

for some xi,X2 and X3(= r). Moreover, e~^2T2g-xir3 ^ ^j^g^^. ^.g^^^ reabsorb it into h 

and write the general point on 5*0(4) as 

g[xi, ...,xe] = e"i^='e"2^2e"^'^^/i[x4,X5,X6]. 

We know the range for the parameters X4, X5, xg. To determine the ranges for the remaining 
parameters, we note that 

p[xi,X2,X3] = e^l^3e^2^^e^•^'^^ 

parametrises the points of 50(4)/50(3) = 5^, so that we need to compute the metric 
induced on 50(4)/50(3) and compare it with the metric of S^. As before, we need to 
compute the left invariant form 

Jp ■= p-^dP. 

However, in general, Jp is not cotangent to P, because P is not a subgroup, so that it will 
have a component tangent to the fiber H. Fortunately, our choice for the product separates 
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P from H orthogonally, and we can obtain the infinitesimal displacement along P by an 
orthogonal projection (that is dropping the terms Tj, i = 1, 2, 3). If we call such a projection 
n, we, thus, have 

dsl = -^Trace(n(Jp)n(Jp)), 

where 

Jp = dxsT4 — sin X3dx2TQ — sin X2 sin x^T^dxi + . . . , 
the dots indicating the terms tangent to H. Thus, we get 

ds^p = dxl + sin^ X3{dxl + sin^ X2dxl), 

which is just the metric of in the usual spherical polar coordinates, 

— {X e R^|Xi = cosxa, X2 — sinxa COSX2, Xs — sinx3sinx2COSXi, X4 — sinx3sinx2sinxi}, 

so that we must choose 

X3e[0,7r], X2e[0,7r], X3e[0,27r]. 
This complete our determination of the ranges. 



The invariant measure is 

d^so{A) — \j d.Qids\j^ deids^p — sin X2 s\v? x^ sin x^dxidx2dx3dx4dxzdxQ. (48) 

Note that there is a second quite interesting way to determine the ranges using the measure. 
To understand it, let us think of a sphere S"^ parametrised with latitudinal and longitudinal 
(polar) coordinates. But suppose we do not know the ranges of parameters. The only fact 
we know is that we can cover our sphere entirely, if any coordinate runs over a whole period 
(A e [0, 27r], (f) e [0, 27r]). The problem is that we can cover the sphere many times (twice 
in our case). But suppose we know the corresponding measure: | sm\\d\d(j). We see that it 
becomes singular at A = and A = tt. These correspond to the points where the parallels 
shrink down (that is, the north and south poles). This implies that when A runs over [0, tt], 
we generate a closed surface which, thus, must cover the sphere (which contains no closed 
surfaces but itself), so that the range of A can be restricted to [0,7r]. 

In the same way, we can look at the invariant metric we have just constructed for 5'0(4). 
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All coordinates xj, I = 1, . . . ,6 have period [0, 2n] because the exp(x/Tf) are functions of 
cos a;/ and sinxj. However, X2, and x^ appear in the measure exactly as for the sphere 
and must be restricted to [0, n]. Thus, we re-obtain the ranges 

xi,X4,,xe e [0,27r], X2,X3,X5 e [0,tt]. (49) 



VI. APPENDIX II (S. CACCIATORI) 

Let us call A the eigenvalues simplex and M the subset of varying angles for the chosen 
example. Finally let us call 5* the subset of M imposed by the Peres-Horodecki condition 



(13) 



- - sin^ (2x5) (Ai - A3) ' - A2 (Ai + A2 + A3 - 1) > 0. 



Thus we need to compute 



J S J M 



(50) 



(51) 



dyU being the measure on the given region. In particular (13) constraints the range of x^ 
only, so that 



ESFi 



1 



sin x^dx^. 



(52) 



Note that (13) is invariant under X5 — > tt — X5 so that we can restrict the region S to 



5*0 = S* n [0, 7r/2] and write 



(53) 



ESF^xs} = / sinx^dx^. 

Jso 

Furthermore, in [0,7r/2] the map = sinx^ is bijective as a map [0,7r/2] — t- [0, 1] so that 



we can write 



where a is the set of solutions of 



ESF{^^} = / dy^ 

J a 



yl{\-yl) (Ai-A3)' + A2A4>0, 



(54) 



(55) 



in [0, 1], A4 = 1 — Ai — A2 — A3. To determine a let us first set z = yf. The solutions of 

(^2_^)(Ai-A3)' + A2A4 = 0, (56) 



are 



1± Wl 



4A2A4 



(Ai - As)^ 



(57) 
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When the discriminant is non positive, that is when (Ai — Aa)^ < 4A2A4, then (55) is satisfied 
for all y and ESF = 1. Otherwise z± satisfy 



< z_ < z, < 1. 



(58) 



In this case 



[z^ - z) (Ai - A3) ' + A2A4 > 0, 



is solved for external values and being z G [0, 1] we have that (55) is satisfied for 



ze [0,z_]U[2;+,l], 



and then (being again a/z = E [0, 1]) 



y, G [0, ^] U [y^, 1] = a. 



(59) 



(60) 



(61) 



Thus 



ESF{x^} = 1 - y/z^ + y/zl =1 j= h 



V2 



V2 



(62) 
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